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ABSTRACT 

Two computer models, the U.S. Environmental Protection Agency's 
SWMM and the U.S. Army Corps of Engineers' STORM, were applied for simulation 
of urban runoff phenomena in a test catchment of the "new town" of Meadowvale 
prior to development and after its completion. A special subroutine for 
the STORM model was developed to simulate the performance of an in-stream 
settling basin-reservoir system for storm water quality control. 

It was predicted that the uncontrolled urbanization of the 
Meadowvale catchment would lead to a deterioration of storm water quality. 
Expected increases in annual pollutant yields (compared with undeveloped 
conditions) would be: BOO from 4 to 26 lb/acre; SS from 3^ to 171 lb/acre; 
nitrogen from 1 to 10 lb/acre, and, orthophosphate from Q.k to I lb/acre. 
Hydrological ly, peak flows might increase by a factor ranging from two 
to five depending on the storm frequency, while the total annual runoff 
would increase from 4.1 inches/annum to 8 inches/annum. The settling 
basin-reservoir system to be constructed at Meadowvale should significantly 
improve runoff characteristics. Predicted yields of pollutants passing 
through the system to the Credit River were: BOD - 1.1 lb/acre; SS - 11 
lb/acre; nitrogen - 7 lb/acre; and orthophosphate - 0.8 lb/acre. Peak 
flow rates would drop, the actual reduction in the flow rate depending 
on the size and temporal pattern of the storm. It was estimated that 
the system designed for Meadowvale would reduce peak flows to pre-urbanized 
levels or less for storms with return periods less than one year. The 
studied system of storm water quality control appears to be very attractive. 
However, future measurements are required to calibrate the applied 
models and verify the accepted assumptions. 



RESUME 



Les auteurs se sont servi de deux modules i nforma t iques , le 
SWHM* de 1'EPA Hp<; F -U. et is STORM du U.S. Army Cui pb of Engineers, 
pour simuler le ru i ssel 1 ement dans un bassin donne d'une ville modele, 
Meadowvale, avant et apres 1 'amenagement. Ms sont mis au point, pour 
!e modele STORM, un sous-programme qui simule l'epuration des eaux 
pluviales dans un bassin de stabilisation. 

II ressort de ces simulations que 1 ' urban! sat ion "sauvage" de 
Meadowvale entrafnerait une forte diminution de la qualite des eaux 
pluviales. Au bout d ' un an, la pollution resultante se traduirait par une 
augmentation de la DBO de 4 a 26 lb/acre, des matieres en suspension (MES) 
de 43 3 171 lb/acre, de I 'azote de 1 3 1 lb/acre, et des orthophospha tes 
de 0,4 a" 1 lb/acre. Du point de vue hydrolog i que, le debit maximal 
pourrait doubler et meme quintupler, selon la frequence des pluies, 
tandis que la quantite d'eau de ru i ssel I ement annuelle passerait de 4,1 
3 8 pouces. L' amenagement du bassin de stabilisation devrait ameliorer 
sensiblement la situation. Apres avoir traverse le systeme, les quant i tes 
prevues de pollutants rejetes dans la riviere Credit seraient les suivantes; 
1,1 lb/acre pour la 0B0, II lb/acre pour les MES, 7 lb/acre pour 1 'azote, 
et 0,8 lb/acre pour les orthophosphates . Les debits maximaux, variables 
suivant 1 ' importance et la duree des pluies, baisseraient aussi. Dans 
le cas de pluies dont le cycle est de moins d ' un an, on estime que le 
systeme concu pour Meadowvale ramenerait les debits maximaux a des valeurs 
egales ou inferieures 5 eel les d'avant 1 ' urbani sat ion . Le systeme 
d'epuration des eaux pluviales etudie semble avantageux. II serait 
toutefois necessaire de prendre d'autres mesures pour etalonner les modeles 
et verifier les hypotheses formulees. 



SWMM (Storm Water Management Model) 
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CONCLUSIONS 

The following conclusions can be drawn from the performed simul- 
ations of the Meadowvale test catchment: 

1) Urbanization of the Meadowvale catchment will lead to a deterioration 
of storm water quality. Expected increases in annual pollutant 
yields as compared with the undeveloped conditions are: 

BOD from k to 26 lb/acre; SS from 3*» to 171 lb/acre; nitrogen 
from 1 to 10 lb/acre; and orthophospha te from 0.4 to I lb/acre. 
Hydrologi cally, peak flows may increase by a factor ranging from 
two to five depending on the storm frequency, while the total annual 
runoff would increase from 4.1 in/annum to 8 in/annum. 

2) The silting basin with the dam to be constructed at Meadowvale 
should significantly improve runoff characteristics. Predicted 
yields of pollutants passing through the system to the Credit 
River are: BOD - 1.1 lb/acre; SS - 1 I lb/acre; nitrogen - 7 
lb/acre; and orthophosphate - 0.8 lb/acre. Peak flow rate drops, 
the actual reduction in the flow rate depending on the size and 
temporal pattern of the storm. It is estimated that the system as 
designed at Meadowvale will reduce peak flows to pre-urbani zed 
levels or less for storms with return periods less than one year. 

3) For the design of storage systems in residential developments such as 
Meadowvale, consideration should be given to the use of two spillways, 
one to restrict flow and allow the use of available storage, and 

the second to cater to large events (i.e., greater than some 
nominal return period, say 25 years). The increased costs of a 
higher dam and higher flood servitudes around the dam should 
not be significant at this scale of project, whereas the benefits 
of peak flow reduction are considerable. It may be possible to 
control storm water flow by designing the system for peak flow 
reduction, the assumption being that if peak flows are reduced to 
pre-urban values, problems of erosion, flooding and pollution 
will be checked at acceptable levels. 
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4) Based on predictions of pollutant concentrations in the reservoir, 
a seasonal algae growth can be expected. This problem should be 
control led. 

5) The expected buildup of sludge will be slow, chiefly due to the 
large area available for sludge storage, averaging 0.05 in/annum. 
Even during construction of the remaining undeveloped land, build- 
up of sludge will only be moderate. The basin will trap flotables 
and has the potential to become anaerobic. A program for regular 
drawdown and cleaning should, therefore, be planned. 
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RECOMMENDATIONS 

1) The storm water quality control system composed of a settlina tank 
followed by a water reservoir, may be recommended as a relatively 
simple and reliable method for storm water pollution abatement in 
new residential developments. 

2) Future monitoring of the Meadowvale test site is required to cali- 
brate the models applied in this study and verify the accepted 
assumptions. Most sensitive parameters are: imperviousness , 
catchment width, infiltration capacity, surface depression, roughness 
and slope (listed in order of decreasing importance). 

3) Many aspects of the suggested storm water control method are still 
not known well enough and require further research. More specif ic- 

al ly: 

(a) The settling behaviour of the storm water solids should be 
studied with particular emphasis on the removals of nutrients 
and heavy metals obtainable under different hydraulic loadings. 

(b) The reservoir quality simulation algorithms should be tested. 
It is probable that a more sophisticated approach (eg. a bio- 
logical model of the nutrient balance) may be necessary if 
acceptable results are to be obtained. 

(c) The significance of snowmelt on the reservoir water quality 
should be studied. 

(d) The microbiological quality of the storm water runoff and the 
effects of prolonged storage on bacterial populations should 
be examined. 



1. INTRODUCTION 

In the past five years there has been a proliferation of 
computer models for the simulation of the urban runoff phenomena. This 
study has been based on two such models: SWMM, developed for the U.S. 
Environmental Protection Agency by Metcalf and Eddy, Inc., the University 
of Florida, and Water Resources Engineers, Inc.; and STORM, developed by 
Water Resources Engineers and currently managed by the U.S. Army Corps 
of Engineers. Both of these models have been used without any initial 
calibration. Therefore, the results presented here have only preliminary 
character and must be verified by future measurements. 

In this report, after presentation of general comments on the 
applied models {Chapter 3), the development of the new subroutine is 
described (Chapter ^). This is followed by the characterization of the 
test site (Chapter 5) and the results of the hydrological (Chapter 6) 
and pollution (Chapter 7) simulations. The new subroutine program 
"MIKE" is given in the appendix to this report. 



2. OBJECTIVES 

The objective of the study reported here was to examine the 
behaviour ot the settling basin-reservoir system used to control quantity 
and quality of storm water from residential areas. For this purpose, the 
computer models and proper data for the test site were used to generate 
prediction of flow hydrographs, pol 1 utographs and annual pollutant 
yields. 

The proposed system at Meadowvale incorporates a siltation basin 
(^0,000 ft /3700 m ) followed by a larger reservoir (10 acresA hectares). 
A more detailed description of the catchment is given in Chapter 5- The 
catchment was simulated for before and after development conditions. 
Particular attention was paid to predicting pollutant yields, the quality 
of the reservoir water and the operational efficiency of the siltation 
basi n. 

Wherever possible, data from Meadowvale were used. Unfortunately 
quality simulations were based on assumed pollutant loading rates. The 
Ontario Ministry of the Environment will be monitoring, over the next 
four to five years, the performance of the Meadowvale system and will 
provide data to check this work. 



3. RUNOFF MODELS 
3- 1 General 

Runoff water quality predictions according to both the applied 
models (SWMM and STORM) are relatively inaccurate. This is due mostly 
to the preliminarily simplistic approach adopted for the logic in these 
models. 

Some discussion of this approach seems essential. 

3.2 SWMM 

The Environmental Protection Agency's Storm Water Management 
Model (SWMM) is probably the most comprehensive urban runoff model 
currently available. It is well documented (Environmental Protection 
Agency, 1971) and is being applied by an increasing number of consultants 
and public authorities. For runoff and flow computations the model 
uses deterministic routines and is known to give acceptable results pro- 
viding no surcharging occurs (MacLaren, 1975). Calibration of sensitive 
parameters is desirable but not essential. 

Quality computations in this model are not as accurate. This 
is partly due to the simplistic approach adopted for modelling runoff 
quality, although the phenomenon itself is complex, stochastic and 
difficult to conceptualise. The salient features of the quality routine 
in SWMM are: 

1) Dust and dirt build up on the subcatchments at a constant 
daily rate which is dependent on gutter length and land 
use. This rate of build up is the most sensitive parameter 
for quality simulation and, in general, the total amount 

of pollutant washed-off during a storm is approximately 
linearly proportional to the number of dry days previous 
to the storm. 

2) Proportionality coefficients and specific pollutant loadings 
of dust and dirt are related to land use and built in as 
default values, making calibration a tedious procedure. 

3) The rate of pollutant washoff explicitly depends on the 
runoff and implicitly on the rainfall intensity, as well 



as on the availability of the pollutant. It is assumed that 
90% of the available pollutant is washed off the catchment 
by a one-hour storm with an intensity of J inch oer hour 
(1.27 cm/hr). No pollutants are generated during runoff. 

A recent paper by Colston and Tafuri (1975) reported the results 
of a study in Durham, North Carolina. Using water quality data from the 
runoff of 66 storms, statistical techniques were used to derive regress- 
ion equations for pollutant concentrations and identify important variables. 
Four variables were investigated: the discharge rate, the time from 
storm start, the time from the previous storm, and the time from the last 
peak. They found that the discharge rate and the time from the storm 
start were significant while the other two variables yielded only a 
small gain in the correlation coefficient and could be ignored. A 
similar finding was also reported in a previous study made in Tulsa, 
Oklahoma (Environmental Protection Agency, 1970). 

Thus, it appears that the SWMM procedure, in which the time 
from the last storm or the number of dry days is a controlling variable, 
may not be precise enough. Possibly the rate of turnover of pollutants on 
a catchment is a much slower process than assumed by SWMM, i.e., a storm 
would only remove a fraction of the pollutants on the catchment (possibly 
the available fraction) but because of erosive, wetting and drying processes 
the available fraction is generated soon after the end of the storm. In 
addition, SWMM does not recognize that for a nonconservative pollutant 
such as BOD, a steady state may exist where the rate of generation 
(from dust and dirt) equals the rate of decay. Conversely, conform 
numbers may in fact increase given suitable environmental conditions. 
It appears that much research is needed to clarify how pollutants are 
deposited, washed-off and transported before a deterministic model can 
be developed. 

Many researchers believe that there is no definable relationship 
for determining the quality of storm water from a specific event and 
that investigations should be directed at predicting aggregate (statistical) 
effects. 



3.3 STORM 

STORM is the planning model of the U.S. Army Corps of Engineers 
(197^)- It permits the comparative analysis of different storm- treatment 
combinations for control of overflows. The program uses a storage-treat- 
ment hierarchy which is shown in Figure 1. This type of priority routing 
is ideally suited to the management of combined sewer overflows; however, 
for separate sewers (the typical arrangement in new residential develop- 
ments), in-stream storage (Figure I) is more practical. Some computational 
facility to simulate an in-stream reservoir is desirable for STORM; such 
a routine was developed and is discussed later in this report. 

The program is designed for use of many years of continuous 
hourly precipitation records, but selected events may also be studied. 
For runoff computations, the program uses the following equation: 

R = C(P - K) • A 

where: R = runoff (cfs), 

C = runoff coefficient, a constant for the catchment, 
P = precipitation in the hour (in), 

K= depression storage at the beginning of the hour (in), 
A * area (acres) . 

Although this appears at first sight to be the "rational formula", 
the precipitation is not based on an intensity-duration frequency curve. 
Depression storage is recovered between storms to a maximum user-specified 
depth. The major criticism of the runoff algorithm is the use of a 
constant runoff coefficient; in practice, the coefficient decreases as 
rainfall accumulates. The quality routine is analagous to that in SWMM 
already discussed; however, two points must be noted: 

1) For SWMM, the initial load of pollutants in the catchment is 
determined by the number of dry days - a user specified 
parameter - but for STORM, buildup is a continuous process 
and in times of no rain may become excessive. The error is 
fortunately damped out. An upper limit for pollution 
buildup is required to solve the problem. 

2) For STORM, the rate of pollutant buildup is an input 
option - a distinct advantage over the SWMM model. 
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SUBROUTINE "MIKE" 

Development of the Subroutine 

This subroutine, developed by M.P.H. Murrey, at the University 



of Toronto as part of this research, is an addition to STORM which permits 
the simulation of an in-stream settling basin-reservoir system. The two 
components are considered independently. A flow chart to illustrate 
the model logic is shown in Figure 2 and a listing of the program source 
and the necessary modifications to STORM are given in the appendix. 

Input data for the subroutine is read in after the "T Group" 
cards described in the STORM User's Manual (197**). Input requirements 
and format are summarized in Table 1. The following sequence of computation 
is made for every hour of the simulation: 

1) If runoff occurs and in-stream modelling is required, the routine 
is called; if not, computation proceeds in Main. 

2) The number of hours since the previous runoff is calculated. 
The reservoir water quality and storage condition is updated 
for each hour using the following relations: 
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where: 



T = time, 

f a -20) 
«! = deoxygenation coefficient = 0.15 x 1 . 0^+7 

(a -20) 

K 2 = reaeration coefficient = 0.075 x 1.0241 v , 
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SAT 



= temperature, C, 

= saturated dissolved oxygen concentration. 
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TABLE I. INPUT REQUIREMENTS AND FORMAT FOR SUBROUTINE MIKE 
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CARD 


FORMAT 


VARIABLE 1 


DESCRIPTION 


1 


15 


IRES 


IRES = NO IN-STREAM RESERVOIR 

>0 IN-STREAM RESERVOIR SIMULATED 


2 


(3F10.1, 215, F6.1) 


PVOL 


VOLUME OF SILTATION POND (ft 3 ) 






ARP1 


AREA OF SILTATION POND (ft*) 






RVOL 


VOLUME OF RESERVOIR FULL (ft 3 ) 






WTP2 


TYPE OF WEIR: 

< INDICATES OGEE 

= INDICATES MORNING GLORY 

>0 INDICATES CONTROLLED DISCHARGE 






NHTS 


NUMBER OF PAIRS OF STAGE AREA DATA POINTS 
TO BE READ IN GROUP THREE 






WDIM 


CHARACTERISTIC DIMENSIONS OF WEIR 
FOR WTP2 < WDIM = LENGTH ft 
FOR WTP2 ■ WDIM - RADIUS ft 
FOR WTP2 > WDIM = FLOW CFS 
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3(F5.1, F8.0) 


HAJ(R) 


STAGE AREA DATA ABOVE 






ARJ(R) 


FULL SUPPLY LEVEL, 3 PAIRS/CARD 
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6(F7.2) 


CONC(K) 


INITIAL VALUES OF PARAMETERS DESCRIBING 
WATER QUALITY 
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For SS and SETL S Equations (1) and (2), 21% reduction in two 
hours is somewhat arbitrarily assumed. Minimum values of 
5 mg/1 and 0.05 mg/l , respectively, are assumed. 
Equations (3), CO and (5) are recognized as the Streeter 
Phelps equations (Velz, 1970). Whether such an approach 
(traditionally used for oxygen budgets in rivers) is valid in a 
reservoir situation is disputable and must await experimental 
evidence. The major cause for this concern is the difference 
in the mixing characteristics of rivers and reservoirs. The 
deoxygena t i on and reaeration coefficients recommended by 
McGauhey (I9&8) are used. The concentration of phosphorus (as 
ortho) and nitrogen are assumed to remain constant between storms 
(Equations (6) and (7)). This is a major simplification since 
factors such as the dissolved oxygen concentration, the presence 
of algae or aquatic vegetation and the chemistry of the sediments 
are ignored. For each hour reservoir storage is incremented 
(Equation (8)) using average monthly potential evaporation rates. 

For the hour of runoff the efficiency of the stilling basin 
(in terms of suspended solids removal) is calculated using 
the equation: 

EFF = 0.9 - 0.2 x 10 x OR 2 , for OR < 1000 U.S. gal/ft 2 /d 
EFF = O . 7e 0.00057(1000 - OR) ^ fm ]Q0Q < QR ^ gal/ft 2 /d 

where: EFF = efficiency, and 
OR = overflow rate. 



The assumed curve is shown in Figure 3- For each pollutant 
the effect of sedimentation is related to the efficiency EFF 
by: 
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These assumed residuals are based on the performance typically 
achieved in primary clarifiers for sewage treatment (Heinke, 
1973, and Metcalf and Eddy, 1972). The flow hydrograph is 
assumed to pass through the siltation basin without attenuation. 

h] For the reservoir, the quality of the water is updated using 
the pollutant inputs from the stilling basin and assuming 
completely mixed conditions. The flow is routed using the 
well known puis technique (Wilson, 1 969) . Three discharge 
options are available, viz. the Ogee type spillway, Morning 
Glory type spillway, and a controlled constant discharge. 

5) A counter variable records, for each month and for each poll- 
utant, the accumulated number of hours in which the pollutant 
concentration has been within a specified range. 

Output from the subroutine includes a frequency analysis of the 
stilling basin efficiency, an estimate of the volume of solids removed 

in the stilling basin and a monthly frequency analysis (for the six 
parameters SS, SETL S, BOD, N, PO,, DO) of the reservoir water quality. 

If detailed event analysis is requested, inflow and outflow hydrographs 
and the concentrations of the six parameters in the outflow are printed. 

The model as developed can and should be further improved. 
Capability for allowing chemical treatment, for including intentional 
drawdown policy, for reading in user options for constants (such as monthly 
temperature, pollutant categories and decay coefficients) and for simulating 
bacteria! contamination would be immediate objectives. It is hoped that 
future work will be done at the University of Toronto (or elsewhere) on 
MIKE and that ultimately the routine will be released for general use as 
part of STORM. 

ft. 2 University of Toronto Experience 

The practical application of STORM requires three temporary data 
sets to store information for the quality report, for the pollutograph 
report and for the erosion report. These temporary data sets are passed 
at the end of the first job step (the STORM computation) to a second 
step which, programmed by the user, must print out the contents of the 
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temporary files. The STORM program does not itself produce this print- 
out. 

Progressing from this stage, rainfall data for the years 1961 - 
197*» (some ten thousand cards) were obtained from the Atmospheric Environment 
Service (AES) of Fisheries and Environment Canada. Many of the data cards 
contained negative rainfalls, alphabetic rainfalls and over-punches and 
the entire set had to be sorted before use. 

The STORM program, as obtained from the Ontario Ministry of the 
Environment, was set up to read United States rainfall cards and before 
the edited AES data could be used, a change of format, presented below, 
was necessary: 

FOR U.S. DATA CARDS 2120 FORMAT (2X, 16, 2^13) 1915 

FOR CANADIAN DATA CARDS 2120 FORMAT (9X, 16, 5X, 1213,/, 20X, 1213) 1915 



5. MEADOWVALE TEST SITE 
5. 1 General 

Meadowvale is a "new town" presently under construction in 
Mississauga, Ontario. The town comprises some 300 acres and the population 
is expected to reach 65,000 after completion in 1985. A simplified plan 
of the development is shown in Figure k. Two dams are planned as features 
of a connected green belt. Construction on the first of these in 
Meadowvale West was to start in 1976. The future lake and its catchment 
is the prototype which has been examined in this study. 

Most of the 2^5~acre catchment will be residential. A portion 
of the catchment is outside of Meadowvale proper (the area to the west 
of Winston Churchill Boulevard). It is assumed that this area will 
ultimately be residential. Prior to development the land was used for 
agriculture; originally, however, it would have been forested. 

The 12 acre dam and stilling basin facility will be used for 
non-body contact activities (boating, fishing etc.). The water quality 
of the dam basin will thus be the major determinant of its recreational 
value. Discharge from the dam is to be channelled to the Credit River, 
approximately two miles distant. 

It is expected that the Meadowvale approach to development 
(i.e., the integration of water-based recreational facilities into 
"completely" planned communities) will dominate planning concepts for 
the next 20 years. 

5.2 Discretization for SWMM and STORM 

The input data requirements for simulation using SWMM are 
complex. The catchment must be discretized and a system of pipes and 
junctions representative of the actual drainage structure must be used 
to connect each subcatchment. Three conditions were modelled using 
SWMM, i.e., undeveloped, developing and fully developed conditions. The 
drainage system designed for Meadowvale will relocate minor watershed 
boundaries, making a comparison of pre and post-development conditions for 
subcatchments difficult. To resolve this problem, pre and post-development 
watershed boundaries for the dam were assumed to coincide with those for 
developed conditions. 



\h 




CATCHMENT BOUNDARY 



FIGURE k. PLAN OF THE MEADOWVALE WEST CATCHMENT 
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For STORM, data requirements are much simpler. The entire 
catchment is lumped and no transport system is used. Conditions before 
and after development (as for SWMM) were simulated and the required 
physical data are shown in Table 2. 

5.3 Rainfall Data 

Rainfall data from the Atmospheric Environment Service (AES) 
station at Mai ton, ten miles from Meadowvale, were used for both SWMM 
and STORM simulations. No areal correction factors were applied and 
snowmelt was not considered. 

For STORM, \k years of hourly rainfall records (1961-1974) 
were used. Since STORM gives a continuous simulation, the effects of 
evaporation on depression storage regeneration must be considered. The 
respective data, also provided by AES, are shown in Table 3. 

For SWMM, continuous recorded rainfall measurements for 1974 
were used to develop ten-minute interval hyetographs. All events (an 
event is defined in this case by 12 hours of no rain prior to the storm) 
with a total rainfall of greater than 0.2" were simulated. Events longer 
than 33 hours, the maximum possible simulation for ten minute hyetographs, 
were contracted. This may give a small overestimate of runoff. Events 
of periods of no rain greater than two hours were contracted to avoid 
difficulties arising from regeneration of excess depression storage 
by infiltration during breaks in rainfall (Orlob, 197*0- 

5.4 Quality Simulation Data 

SWMM and STORM use analogous procedures for determining runoff 
quality. The required data are: 

1) the daily rate of dust and dirt accumulation in lb/100 ft 
of curb for each land use; 

2) the pounds of specific pollutant per 100 lb of dust and 
dirt (STORM has the facility to model nonurban areas, in 
which case loadings are reported directly In pounds of 
pollutant per acre); and, 

3) catchment parameters such as the length of gutters, the 
number of catch basins, the frequency of streetsweep ing , 
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TABLE 2. PHYSICAL CATCHMENT DATA FOR STORM SIMULATIONS 



1. UNDEVELOPED CONDITIONS 

NON-URBAN WATERSHED DATA 

LAND AREA - 241 ACRES (97 HECTARES) 
RUNOFF COEFFICIENT = 0.25 
DEPRESSION STORAGE = 0.1 IN (2.5mm) 

URBAN WATERSHED DATA 

LAND AREA = 1 ACRE 

LAND USE = SINGLE FAMILY RESIDENTIAL 
RUNOFF COEFFICIENT 

(PERVIOUS) = 0.15 (75% OF AREA) 

(IMPERVIOUS) - 0.9 

DEPRESSION STORAGE « 0.05 IN (1.3mm) 

GUTTER LENGTH = ? FT/ ACRE (1.5m/Ha) 



URBAN WATERSHED DATA 

LAND AREA = 

RUNOFF COEFFICIENT 

(PERV AREA) = 

(IMPERV AREA) = 

DEPRESSION STORAGE ■ 



2. DEVELOPED CONDITIONS 

244 ACRES (98 HECTARES) 

0.15 

0.9 

0.05 IN (1.3mm) 



LAND USE 



% OF 
LAND AREA 



PERCENT IMPERV 



GUTTER LENGTH 
FT/A M/Ha 



SINGLE 
MULTPL 
COMMCL 



65.2 
22.9 
11.9 



30.0 
65.0 
85.0 



300.0 226.0 
300.0 226.0 
400.0 301.2 
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TABLE 3. MONTHLY AVERAGE EVAPORATION RATES AND TEMPERATURES FOR 
MALTON WEATHER STATION (AES) 



MONTH 


MEAN TEMP 


MEAN 


EVAP 




°C 


inches 


mm 


JANUARY 











FEBRUARY 











MARCH 


6.3 


0.3 


8 


APRIL 


6.4 


3.6 


91 


MAY 


12.0 


4.8 


122 


JUNE 


18.2 


5.7 


145 


JULY 


20.7 


6.9 


175 


AUGUST 


20.0 


4.8 


122 , 


SEPTEMBER 


15.7 


3.6 


91 


OCTOBER 


9.8 


1.8 


46 


NOVEMBER 


4.5 


1.2 


30 


DECEMBER 
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the efficiency of streetsweeping and (for SWMM) the number 
of dry days prior to the event. 

Data on loading rates of dust and dirt and on the component 
fractions of the dust and dirt were not available for Canadian conditions; 
hence, the values shown in Table h, which are based on the APWA study in 
Chicago (American Public Works Association, 1 969) , were assumed. Work 
is presently being done under the aegis of the Canada-Ontario Agreement 
to establish this information for Canada and should be completed shortly. 

Quantitative values for the various catchment parameters were 
obtained from a variety of sources, including the Public Works Department 
of the City of Mississauga, Project Planning Associates Limited (the 
project consultants) and a personal inspection of the site. 

The regression equations developed by Colston and Tafuri (1975) 
were used to give a second estimate of runoff quality. Flow hydrographs 
obtained from SWMM provided the necessary input data (i.e. flow rate and 
time from start of storm). 

5- 5 Reservoi r Data 

Input data for the reservoir simulation using subroutine MIKE 
are presented in Table 5- The initial quality of dam water was arbitrarily 
chosen. For a long-term simulation, such as was done for Meadowvale, 
the choice is unimportant. 



19 



TABLE k. ASSUMED POLLUTANT LOADING RATES IN SWMM AND STORM 



O 



LANDUSE 


D&D RATE 
LB /lOO'/D 


LBS POLLUTANT / 100 LB OF D&D a 


COLIFORM 
MPN/gm OF D&D 


SUSP S c 


SETL S C 


B0D C 


COD 


N C 


po 4 c 


GREASE 


SINGLE 

FAMILY 


0.7 


11.1 


1.1 


0.5 


4.0 


0.048 


0.005 


0.1 


1.3xl0 6 


RES. 




















MULTPL 
FAMILY 


2.3 


8.0 


0.8 


0.36 


4.0 


0.061 


0.005 


0.1 


2.7xl0 6 


RES. 




















COMMER- 


3.3 


17.0 


1.7 


0.77 


3.9 


0.041 


0.007 


0.1 


1.7xl0 6 


CIAL 




















INDUST- 


4.6 


6.7 


0.7 


0.3 


4.0 


0.043 


0.003 


0.1 


l.OxlO 6 


RIAL 




















OPEN- 


1.5 


11.1 


1.1 


0.5 


2.0 


0.005 


0.001 


0.1 





PARK 




















NON . 
URBAN 


N/A 


0.023 


0.02 


0.03 


- 


0.01 


0.003 


- 


- 





















a MOST VALUES ARE BASED ON 1969 APWA REPORT 

b THIS LANDUSE IS AN OPTION IN STORM ONLY, UNITS ARE LBS/ ACRE/DAY. 

c POLLUTANTS SIMULATED IN STORM. 

CONVERSION FACTORS: 1LB/100 FT ■ 15gm/m, 1LB/100LBS * lOmg/gm, 1LB/ACRE = 1. lKg/HECTARE 



TABLE 5. INPUT DATA FOR RESERVOIR SIMULATION USING SUBROUTINE MIKE 



1. SILTATION BASIN 

VOLUME AT WHIR LEVEL = 474200 FT 3 (13427m 3 ) 
AREA AT WEIR LEVEL = 39,520 FT 2 (3671m 2 ) 

2. DAM ** ** 

VOLUME AT F.S.L. (928.2') = 4.85xl0 6 FT 3 (1.37xl0 5 m 3 ) 

HEIGHT ABOVE WEIR AREA 

2 



FT 


1.8 
3.8 



m 


0.54 
1.15 



455200 
481920 
512480 



m 

42289 
44772 
47610 



3. SPILLWAY 

DROP-INLET TYPE 

CHARACTERISTIC DIMENSION = RADIUS = 8' (2.4m) 



■ 4. INITIAL DUALITY OF DAM WATER 


CONC (1) = S.S. 


= 20mq/l 


CONC (2) = SETL. S. 


- 5mg/1 


CONC (3) = BOD 


= 10mq/l 


CONC (4) = N 


= lmq/1 


CONC (5) = P0 A 


= 0. lmq/1 


CONC (6) = DISSOLVED 


OXYGEN = 8mq/l 
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6. HYDROLOG I CAL SIMULATION 

6. 1 Genera I 

The SWMM model was run for urbanized conditions for all storms 
of 197^ with total rainfalls greater than 0.2" (5 mm). The statistics 
for each storm are listed in Table 6. From this series, three storms 
were chosen for detailed analysis, i.e., to simulate undeveloped conditions 
(allowing comparative analysis of urbanization effects), to study the 
advantages of storage and to investigate erosion. These storms were: 

1) Storm of 4th March - A long (13-hr) frontal storm (1.2"/30 min). 

2) Storm of 1 6th June - The largest (10-hr) storm of the year 

(2.25'757 mm). 

3) Storm of 3rd July - A short convective thunderstorm (1. 3"/ 

33 mm in 3 hr) with the largest ten-minute 
intensity (3" per hour/76 mm per hour). 

The STORM model was run for all storms from 1961 - 197**, inclusive, for 
both developed and undeveloped conditions. 

Comparison of predicted hydrographs using SWMM and STORM are 
presented in Figures 5~10. For urbanized (developed) conditions, two 
features were noted: 

1) Peak flows predicted by SWMM were approximately k5% greater 
than those predicted by STORM. This result was not surprising 
since STORM, by using hourly rainfall data, must consistently 
underestimate peaks from shorter duration effects. The model 
however, was not designed for accurate event prediction. 

2) The total runoff predicted by STORM for each storm was 
approximately 10% greater than that predicted by SWMM. This 
very small difference is a calibration problem requiring 
field data for its solution. Balancing of the SWMM and 
STORM results would magnify the peak discrepancy discussed 
above. 

For nonurbanized (undeveloped) conditions, peak flows as well as 
total runoffs predicted by STORM were significantly greater than those 
predicted by SWMM. A detailed analysis of this is given in Table 7- 
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TABLE 6. STATISTICS FOR 197^ STORMS SIMULATED USING SWMM 



STORM OF 


TOTAL 
IN 


RAINFALL 
mm 


MAXIMUM 10 MINUTE 
IN/HR 


INTENSITY 
mm/HR 


DURATION 
HOURS 


JAN 22 \ 


0.21 


5.3 


0.18 


4.5 


7 


JAN 26 


0.64 


16.2 


0.78 


19.8 


4 


MAR 4 


1.20 


30.4 


0.48 


12.2 


13 


APR 1 


0.48 


12.2 


0.84 


21.3 


2 


APR 3 


1.08 


27.4 


0.6 


15.2 


6 


APR 12 


0.28 


7.1 


0.36 


9.1 


5 


APR 14 


0.48 


12.2 


1.14 


29.0 


3 


APR 22 


0.51 


13.0 


0.36 


9.1 


6 


MAY 5 


0.26 


6.6 


0.30 


7.6 


3 


MAY 8 


0.78 


2.0 


0.18 


4.6 


14 


MAY 15 


0.21 


5.3 


0.36 


9.1 


3 


MAY 16 


2.25 


57.2 


2.46 


6.2 


10 


MAY 28 


0.43 


10.9 


0.54 


1.4 


7 


JUN 13 


0.34 


8.6 


0.42 


1.1 


4 


JUN 15 


0.90 


22.9 


1.98 


50.3 


9 


JUN 19 


1.20 


30.5 


0.96 


24.4 


13 


JUN 21 


0.65 


16.5 


1.5 


38.1 


2 


JUL 3 


1.30 


33.0 


3.03 


76.9 


3 


JUL 29 


0.23 


58.4 


0.36 


9.1 


5 


JUL 30 


0.52 


13.2 


0.90 


2.3 


5 i 


SEP 2 


0.44 


11.2 


0.18 


4.6 


11 


SEP 12 


0.26 


6.6 


0.3 


7.6 


4 


SEP 28 


0.76 


19.3 


0.24 


6.1 


14 


OCT 14 


0.32 


8.1 


0.30 


7.6 


5 


NOV 5 


0.47 


11.9 


0.42 


10.7 


6 


NOV 12 


0.44 


11.1 


0.12 


3.0 


12 


NOV 13 


0.28 


7.3 


0.24 


6.1 


4 


NOV 20 


1.20 


30.5 


0.36 


9.1 
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TOTAL SIMULATED RAINFALL = 18.1" (460mm) 

TOTAL RECORDED RAINFALL = 25.1" (638mm) (FOR 1974) 
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TABLE 7. COMPARISON OF PREDICTED HYDR0GRAPH5 FOR UNDEVELOPED CONDITIONS 
USING SWMM AND STORM 



STORM 


■ 
MARCH 4 


■ — - —■ — ■ — 

MAY 16 


JULY 3 


DURATION (HOURS) 


13 


10 


3 


RAINFALL (IN) 


1.19 


2.25 


1.3 


SWMM: 








RUNOFF (IN) 


0.11 


0.25 


0.21 


PEAKFLOW (CFS) 


7.0 


33.0 


43.3 


STORM: 








RUNOFF (IN) 


0.29 


0.53 


0.30 


PEAKFLOW (CFS) 


12.8 


50.1 


59.3 


RUNOFF 

storm 

RUNOFF 


2.8 


2.1 


1.5 


^wmm 








PEAK 

storm 


1 .8 


1 .5 


1 4 


PEAK 








swmm 
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Not only did STORM predict a greater runoff yield, but the 
prediction was different for each storm, the maximum difference occurring 
with the longest storm (March Ath) . An explanation of this apparent 
anomaly is to be found in the way the two models consider losses. 

STORM assumes a depression storage and any rainfall in excess of 
this storage is available for runoff, but not necessarily converted into 
runoff dependent on the value of the runoff coefficient. SWMM computes 
losses due to depression and infiltration explicitly. Thus, if the 
assumed infiltration rates used in SWMM are too large, an underestimate 
of runoff approximately proportional to the storm length as observed 
would be expected. In addition, the runoff coefficient assumed for STORM 
may be too large. Again this is a calibration problem requiring 
adjustment of model parameters (e.g. infiltration rates, percent imper- 
viousness for SWMM, and runoff coefficients and depression storage for 
STORM) to fit measured flows. Once the runoff predictions agreed, the 
peak ratio shown above would likely become fractional, i.e., STORM 
would underestimate peaks. 

The predictions were considered acceptable for a comparison of 
the hydrological effects of urbanization and storage. 

6. 2 Hydrological Effects of Urbanization 

Comparative studies of urban versus rural drainage basins 
indicate that as the relative magnitude of flood peaks increases, the 
ratio of urban peak flow rate to rural peak rate declines, the effect of 
urbanization being more pronounced for the more frequent occurrences. 
For large storms, infiltration losses are a small proportion of the total 
rainfall and the catchment behaves as if it were impervious. 

At Meadowvale the effects of urbanization as predicted by 
STORM are summarized in Table 8. The annual runoff is seen to increase 
by a factor of 2.15. For specific events, the ratio of nonurban peak 
to urban peak (for 15 large events) was calculated to be constant at 
1.95, the same as the ratio of the runoff coefficients. This is not 
surprising since STORM uses a constant runoff coefficient; however, it 
means that one cannot derive any relationship between storm frequency 
and the peak flow increase due to urbanization using STORM results. 



31 



TABLE 8. STORM PREDICTIONS PRE AND POST-DEVELOPMENT 

Average Annual Statistics for 14 Years of Record 
For the Period Beginning Apr. 9, 1961 and Ending Dec. 12, 1974 

Non-Urbanized Conditions 

Number of Events = 69. 
Number of Overflows = 69-0 

Total Precipitation on Watershed 21.67 inches 

Total Runoff from Watershed 4.09 inches 

Overflow to Receiving Water 4.10 inches 

Initial Overflow to Receiving Water 2.35 inches 

Fraction of Rainfall = 0.19 

Fraction of Rainfall = 0.19, of Runoff = 1.00 

Fraction of Rainfall = 0.11, of Runoff - 0.57 

Average Annual Statistics for 14 Years of Record 
For the Period Beginning Apr. 9, 1961 and Ending Dec. 12, 1974 

Urbanized Conditions 



Number of Events = 69. 1 

Number of Overflows = 69.1 

Total Precipitation on Watershed 21.67 inches 
Total Runoff from Watershed 8.85 inches 

Overflow to Receiving Water 8.85 inches 

Initial Overflow to Receiving Water 5-47 inches 

Fraction of Rainfall = 0.41 

Fraction of Rainfall = 0.41, of Runoff = 1.00 

Fraction of Rainfall = 0.25, of Runoff = 0.62 



The reason for the difference in the peak flow rate ratio and 
the runoff ratio is the depression storage. The depression storage 
assumed (0.1" for undeveloped and 0.05" for developed conditions) determines 
a critical value for total storm rainfall below which no runoff occurs. 
This implies that for some storms between (0.05" and 0.1") there is 
effectively an infinite increase in runoff. 

Using SWMM, the effect of urbanization on specific events is 
illustrated in Figures 11-13. For the three storms, the following 
information was derived. 



STORM 
OF 


TOTAL 
RAINFALL 
in 


15 MIN. 

RAINFALL 
i n/hr 


R 


. P . * 


BEFORE 


URB. 


AFTER 


URB. 


R A 


P A 


/ 




RUNOFF 

in 


PEAK 
cfs 


RUNOFF 
in 


PEAK 
cfs 


R B 


P B 


MHm 


1.2 


0.44 




1 


0.11 


7 


0.50 


34.5 


4.3 


4.9 


16/5/74 


2.25 


1.96 




1 


0.25 


33 


0.94 


135.1 


3.8 


4.1 


3/7/74 


1.3 


2.56 




3 


0.21 


45 


0.55 


193.7 


2.6 


3.9 



Return periods based on maximum 15 minute rainfall intensity and the 
IDF curve for Mai ton 



One weakness of such a comparison is that the storms have different 
temporal patterns. Ideally, to study the effects of urbanization on storms 
of different return periods, linearly proportional storms should be used 
(i.e., storms should be of the same length and pattern, the only difference 
being in the rainfall intensities). 

Generally, these results followed the expected pattern; the effects 
of urbanization were most pronounced for the March **th storm. A study 
of the literature indicates that the predicted increases in peak flow 
may be excessive. Rantz (Shubinski, 1974), using regression analysis to 
derive flood frequency equations, developed a factor to account for 
urbanization. For catchments similar to Meadowvale (i.e., k0% impervious 
and completed sewered), he predicted an increase of 2.5 in the peak flow 
for a one-year storm, compared to the 3-8 predicted by the SWMM for 
the July 3rd storm. In a survey on the effects of urbanization on run- 
off, Rao et al (1972) found increases in peak flow ranging from 1.5 to 
3.0. 
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This meant that either the results for urbanized conditions 
were too large or that the results for undeveloped conditions were too 
small. It was concluded that the SWMM predictions for undeveloped 
conditions were too low for two reasons: 

a) SWMM is designed for urban runoff computations and has 
not been widely tested on undeveloped catchments. 

b) The STORM prediction of runoff using a relatively low 
coefficient of 0.25 was between 50% and 100% greater than 
the SWMM result. 

6.3 Storage Effects 

Net evaporation from the Meadowvale dam (l0"/year) represents 
approximately O.h" of runoff from the catchment or a reduction of 5% in 
the annual yield of 7.08" (corrected average annual runoff predicted by 
STORM). Therefore, the storage system at Meadowvale will have little 
effect on runoff yield. The influence on peak flow is, however, much 
greater. Figures 11-13 show how, for specific events, the incoming 
hydrograph is attenuated and discharged over a much longer period. 

At the outset of this work, it was hoped that the Meadowvale 
system would give a reduction of peak flow to pre-urbanized levels for 
storms up to some nominal return period (10 years). Apparently, this 
level of reduction is only obtained for relatively small storms, e.g., 
Ath of March. Considering the physical size of the dam, with an area of 
± 10 acres and available storage (between the Morning Glory and the 
Emergency Spillway levels) equivalent to 1.5" of runoff, there is obviously 
scope for improvement of the discharge characteristics. 

The basic reason for the poor performance (in terms of peak 
reduction) is the eight-foot radius drop-inlet spillway. Designed using 
the rational method, to discharge "Hurricane Hazel" with a maximum head 
on the crest of 2.1', the spillway does not throttle low flows and 
prevents efficient use of storage. For the May 16th storm, for instance, 
the predicted maximum level reached was only 7". While there is certainly 
a need to safely pass the possible maximum flood, it is noted that an 
emergency spillway, provided at a level three feet above the drop-inlet 
spillway, is adequate for this purpose. 
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To investigate possible improvements in the spillway design, 
three alternatives were examined: drop-inlet spillways with radii of 
two feet and four feet, and a controlled outlet (e.g., pumping) of 8 
cfs. For the three large storms of 197*+, the discharge hydrographs are 
shown in Figures 14-16. In addition, the largest storm for each year 
from 1961-1974, inclusive, was studied using the STORM predictions of 
flow hydrographs. The flows were routed and the maximum discharge and 
water level reached for each storm was calculated. The computation was 
not performed for the constant discharge cases. These results are 
examined in Table 9- 

The following observations are relevant: 

1) Both models first calculate the pollutant mass washed off 
during a time interval and this mass is then divided by 
the runoff volume to obtain concentration of the pollutant. 
In many cases, concentration of the pollutant is propor- 
tional to the flow rate, but may also decrease with the 
increasing flow rate. 

2) The use of a spillway as a flow restriction results in a 
considerable reduction in peak flows for a relatively 
small change in storage utilization. 

3) For a controlled constant discharge of 8 cfs (which incidently 
would allow 3 feet of drawdown in *+5 hours), the maximum 

head reached for the May 16th storm was 13", i.e., by 
increasing the storage utilization by 6", the peak flow is 
reduced from 65 to 8 cfs. 

4) The system as designed at Meadowvale, if subjected to the 
rainfall history for the period 1961-197^, would record 

a maximum water level above the spillway of one foot, 
corresponding to a peak discharge of 1 86 cfs. With a 
four-foot radius drop- inlet spillway, these values change to 
1.3 feet and 1^0 cfs, respectively. 

5) Generally, the relative reduction in peak flows due to the 
dam (measured as ratio of peak flow without the dam to peak 
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TABLE 9. ANALYSIS OF THE EFFECTS OF VARIOUS DISCHARGE OPTIONS ON OUTFLOW CHARACTERISTICS 





STORM 


RANK 


PEAK 9 PEAK 3 

BEFORE AFTER 

URB URB 

P1{CFS) P2{CFS) 


P2 

pt 


STORAGE EFFECTS 






8' 


RAD SPILLWAY 




4' RAD SPILLWAY 




2' 


RAD SPILLWAY 






PEAK 
P3(CFS) 


P3 
PT 


P3 C 
FT 


MAX 

H" 


PEAK 
P3(CFS) 


P3 P3 C 

PT PT 


MAX 

H" 


PEAK 
P3(CFS) 


P3 

PT 


P3 C 
PT 


MAX 
H" 




31/7/70 


1 


107 209 


1.95 


185 


.89 


1.74 


12 


139 


.66 1.28 


16 


64 


.31 


.60 


21 




22/8/68 


2 


95 166 


1.95 


1*5 


.89 


1.74 


11 


125 


.67 1.30 


15 


63 


.33 


.64 


19 




24/7/69 


3 


62 122 


1.96 


109 


.89 


1.74 


8 


85 


.70 1.37 


11 


51 


.42 


.82 


14 




16/5/74 


4 


50 97 


1.94 


86 


.89 


1.74 


7 


68 


.70 1.35 


9h 


45 


.46 


.90 


m 




12/7/62 b 


5 


44 89 


2.02 


63 


.71 


1.43 


6 


41 


.46 0.93 


7 


24 


.27 


.54 


8 




2/4/67 


6 


41 80 


1.95 


60 


.75 


1.46 


$k 


41 


.51 0.99 


7 


24 


.30 


.58 


8 




2/8/65 


7 


37 72 


1.95 


54 


.75 


1.46 


5 


37 


.51 0.99 


6«s 


22 


.31 


.60 


8 


-C- 


3/9/66 


8 


36 70 


1.94 


54 


.77 


1.50 


5 


37 


.53 1.03 


6H 


22 


.31 


.60 


7^ 


ro 


28/8/64 


9 


35 69 


1.97 


51 


.74 


1.46 


5 


35 


.51 1.00 


6 


22 


.32 


.63 


75* 




Z6/7/71 b 


10 


34 73 


2.14 


50 


.68 


1.45 


5 


32 


.44 0.94 


54 


18 


.25 


.54 


6is 




20/7/63 


11 


33 64 


1.94 


42 


.66 


1.28 


4 


1 27 


.42 0.81 


5 


18 


.28 


.54 


6*5 




2/10/73 


12 


30 58 


1.93 


42 


.72 


1.39 


4% 


28 


.48 0.93 


5*5 


16 


.28 


.54 


6 




24/9/61 b 


13 


26 56 


2.15 


36 


.64 


1.38 


4 


22 


.39 0.84 


4% 


12 


.21 


.45 


5 




27/6/72 b 


14 


16 49 


3.06 


29 


.59 


1.80 


3% 


18 


.37 1.13 


4 


10 


.20 


.61 


4% 



NOTES: 



a THESE INFLOW HYOROGRAPHS ARE BASED ON PREDICTIONS BY STORM AND ARE LOW. 

b THESE STORMS EXHIBIT LARGE P2/P1 RATIOS BECAUSE PEAK FLOWS OCCURRED BEFORE DEPRESSION STORAGE WAS FULL. 

c THIS RATIO IS INDICATIVE OF THE OVERALL EFFECTS OF URBANIZATION AND STORAGE ON STORM PEAKS. 



flow with the dam) is greater for the smaller, more frequent 
storms. 

6) For the design of storage facilities such as the Meadowvale 
dam, consideration should be given to the use of two 
spillways, a lower spillway designed to restrict flows to 
pre-urbani zed levels for a nominal return period, and an 
upper spillway to function during low frequency (emergency) 
floods. For large dams this procedure is obviously 
uneconomic, requiring expensive flood servitudes and a 
higher dam, but for smaller systems these are not usually 
problems. For Meadowvale, a four-foot radius drop-inlet 
spillway in conjunction with an emergency spillway gives 
considerable improvement in flood characteristics. 

7) The results in Table 9 are based on STORM predictions which, 
as has been discussed, yielded overestimates of total runoff 
and underestimates of peak flows. Although these factors 
tend to cancel each other, the values given here are probably 
correct only to within 20% error. 
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7. POLLUTION SIMULATION 
7. 1 Genera 1 

As previously discussed, the available storm water simulation 
models are not yet developed enough to give accurate predictions of the 
quality of storm water even if calibration data are available. Given 
the complexity of the processes involved and the difficulty of obtaining 
useful data, the full development of such a predictive capability in 
the near future is unlikely and the present simplistic models will 
continue to be an important tool for the evaluation of storm water 
pollution problems. 

Modelling of runoff quality for undeveloped catchments using 
SWMM is extremely difficult as an arbitrary gutter length (pollutant 
source) has to be chosen (in this case a constant 300 feet for each sub- 
catchment). Presented in Figures 17 and 18 are predicted pol lutographs 
using SWMM for both urbanized and undeveloped conditions (units are 
given both as concentration (mg/1) and mass emissions (lb/min); generally, 
for storm water studies, mass flow units which include the effects of 
a highly variable flow are preferred). The mass pol lutographs for unde- 
veloped conditions have been adjusted up by a factor of two to account for 
the suspected discrepancy in the runoff. Erosion was considered only in the 
prediction for undeveloped conditions as it did not seem important for 
a developed catchment. 

The following observations were made: 

1) A severe degradation in the quality of the runoff due to 
urbanization is predicted. BOD, N and P0, yields will 
increase approximately thirty times while SS yields will 
increase up to fifteen times. The difference between the 
BOD and SS increases is due to erosion, which was considered 
for the undeveloped condition. Basically, the conclusion 

Is that developed conditions without erosion considered 
produce more SS than undeveloped conditions with erosion 
considered. 

2) A comparison with typical measured pollutant concentrations 
in storm water (Table 10), indicates that for developed 
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SWMM POLLUTOGRAPHS FOR THE STORM OF MAY 16, 1974, BEFORE AND 
AFTER DEVELOPMENT. B includes erosion; A does not include 
erosion. 
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TABLE 10. COMPARISON OF THE QUALITY OF STORM SEWER DISCHARGES. 

(DATA ABSTRACTED FROM ENVIRONMENTAL PROTECTION AGENCY, 197*0 
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CITY 


BOD 
mg/1 
AVE RANGE 


SS 
mg/1 
AVE RANGE 


CO LI FORMS* 
MPN/100 ml 
AVE RANGE 


TOT N 

mg/1 as N 

AVE 


TOT P 
mg/1 as P 
AVE 


ANN ARBOR 

MICHIGAN 

1965 


28 11-62 


2080 650-11900 


~ — 


3.5 


1.7 


DES MOINES 

IOWA 

1969 


36 12-100 


505 95-1050 


"" — 


2.2 


0.87 


DOR HAM 

N. CAROLINA 

1968 


31 2-232 


™ — 


3 x 10 5 3 x I0 3 - 2 x 10 


~ 


™ 


SACREMENTO 
CALIFORNIA 
1968-69 


106 2V283 


71 3-211 


S 4 7 

8 x 10 3 2 x 10 - 1 x I0' 


~ 


— 


TULSA 

OKLAHOMA 

1968-69 


11 1-39 


2^7 8^-2052 


6 x 10 5 1 x I0 5 - 3 x 10 


0.3 - 1.5** 


0.2 - 1.2** 


WASHINGTON 

D.C. 

1969 


19 3-90 


1697 130-11280 


1 x 10 5 1 x 10 3 - 5 x 10 8 


2.1 


M 


WINDSOR 
ONTARIO 
1972 


12 2-52 


305 2-M22 


2 x 10 6 1 x I0 3 - 2 x 10 7 


1.6 


0.93** 


MEADOWVALE 
PREDICTED 
(SWMM) 
DEVELOPED 


- 1-50 


5-500 


1 x 10 3 - 1 x 10 6 







only fecal 

nitrogen as organic (kjeldahl), phosphorus as orthophosphate 



conditions the predictions are of the correct order of 
magnitude. The poi lutographs for undeveloped conditions 
appear, however, to be too low, possibly due to an incorrect 
choice of gutter length. 

3) The coliform predictions exhibited computational instability 
(Figure 18). This may have been due to insufficient 
machine storage for the very large number being modelled. 

h) The suspended solids prediction for undeveloped conditions 
was not sensitive to the washoff of dust and dirt {i.e. 
no first flush effect) because erosion was the major component 
of the suspended solids. 

For flooding considerations, the design event is invariably a 
low-frequency storm (e.g., a ten-year storm for urban drainage design), 
but for quality considerations the average monthly storm may be more 
important In terms of the mass of pollutants contributed than the ten-year 
storm. Therefore, for the management of storm water quality, the concern 
seems to be with aggregate pollutional loadings and not with the effects 
of specific events. 

7.2 Pollutant Yield Predictions 



The University of Toronto version of STORM (subroutine MIKE) 
was used to predict annual pollutant yields from the Meadowvale catchment 
to the Credit River for three conditions, i.e., an undeveloped catchment, 
an urbanized catchment, and an urbanized catchment with runoff routed through 
a settling basin-reservoir system. These predictions are summarized in 
Table 11. If the settling basin-dam system for Meadowvale behaves as 
assumed, the yields of SS, SETL S and BOD will be reduced below levels 
existing before urbanization. The system does not have much effect on 
soluble conservative pollutants; nitrogen and orthophospate are reduced 
only by 25% and 16%, respectively. However, it is not expected that either 
of the two models studied will give good estimates for predevelopment 
condi t i ons . 

Shown in Table 12 are pollutant yields reported by other 
researchers; the STORM predictions are seen to be of the correct magnitude. 
Comparing the average pollutant load discharge from Meadowvale with that 
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TABLE II. PREDICTED POLLUTANT YIELDS FOR MEADOWVALE USING STORM 

Average Annual Statistics for 14 Years of Record 
For the Period Beginning April 9, 1961, and Ending December 12, 1974 

Before Urbanization SUSPS SETL S BOD N PO, 

Total Pounds Washoff From Watershed 

Total Pounds Overflow to Receiving Water 

Concentration of Pollutants in Overflow to Receiving Water (mg/l) 

Fraction of Total Washoff Overflowing to Receiving Water 

Fraction of Total Washoff Initially Overflowing to Receiving Water 

After Urbanization SUSPS SETL S BOD N PO 



8131 


699 


1055 


232 


87 


8128 


706 


1060 


247 


106 


36.16 


l,\k 


4.72 


1 .10 


0.47 


1.00 


1.01 


1 .00 


1 .06 


1.22 


0.68 


0.68 


0.68 


0.72 


0.82 



Total Pounds Washoff from Watershed 41409 3165 6239 2284 211 

Total Pounds Overflow to Receiving Water 41404 3169 6241 2289 231 

Concentration of Pollutants in Overflow to Receiving Water (mg/l) 84.61 6.48 12.75 4.68 0.47 

Fraction of Total Washoff Overflowing to Receiving Water 1.00 1.00 1.00 1.00 1.09 

Fraction of Total Washoff Initially Overflowing to Receiving Water 0.70 0.68 0.73 0.71 0.77 

Total Pollutant Yields Overflowing from Dam (lb) SUSPS SETL S BOD U PO^ 

2658.9 67.9 299-8 I678.7 188. I 

Total Pollutant Yields Overflowing from Settling Basin (lb) SUSPS SETL S BOD N^ P0, 

7122.7 720.6 3616. 7 1718.7 193.0 



TABLE 12. COMPARISON OF REPORTED POLLUTANT YIELDS FOR DIFFERENT 
CITIES WITH PREDICTED YIELDS FOR MEADOWVALE 



Locat ion 


LOADINGS 
BOD 


IN POUNDS/ACRE/YEAR* (kg/ha/yr) 
SS Total N 


Total P 


OXHEY 
ENGLAND 


4 
(4.5) 


103 
(115) 


- 


~ 


CINCINNATI 
OHIO 


33 
(37) 


730 
(818) 


8.9 
(9.9) 


2.5 

(2.8) 


ANN ARBOR 
MICHIGAN 


31 

(35) 


1010 
(1132) 


1.9 
(2.1) 


2.8 
(3.1) 


DURHAM 

N. CAROLINA 


34 
(94) 


- 


- 


3.4 
(3.8) 



WINDSOR 
ONTARIO 



MEADOWVALE 
PREDICTION. 
DEVELOPED 
CONDITIONS 
USING STORM 



22 
(25) 



23 
(26) 



440 
(493) 



154 
(172) 



2.85 
(3.19) 



8.4 
(9.4) 



7 - 96-' 

(3.9) 



0.78** 
(0.87) 



• These data abstracted from references 
"-• Or triphosphate only 
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from a typical secondary treatment plant serving the same population, the 
following resutls were obtained: 

Estimated population in Meadowvale catchment at 20/acre = 4800. 

- At BOD equivalent = 0.17 Ib/person/day , BOD load = 300,000 lb/year. 

- At SS equivalent = 0.2 Ib/person/day, SS load = 350,000 lb/year. 

- Assuming a 901 treatment efficiency, BOD yield = 30,000 lb/year, 
compared to 5600 lb from storm water, and SS yield = 35,000 lb/year 
compared to 37,000 lb from storm water. 

In terms of BOD, Meadowvale storm water is an insignificant nonpoint 
pollutant source; however, yields of suspended solids and nutrients can 
clearly not be ignored. Considering the quality of the Credit River 
(data for the years 1 965 _ 7 1 are given in Table 13), it is apparent that 
continued urbanization of the basin, with the associated increases in 
pollutant inputs (particularly nutrients and suspended solids), will 
damage the river and reduce its worth. The above analysis has ignored 
contaminants such as chlorides, asbestos, heavy metals (predominantly 
lead) and PCB's, which are known to occur in urban storm water. 

Pollutant yields for 1974 (for urbanized conditions) were 
also calculated using the SWMM predictions and the regression equations 
of Colston and Tafuri (1975). This information is given in Table ill. 
The predictions using the three methods are conspicuously different and 
some explanation is required if the previous discussion is to have any 
credibi 1 i ty. 

Firstly, the regression equations developed by Colston and Tafuri 
are valid for a specific catchment in North Carolina. Although this 
catchment has similar land uses to Meadowvale, the drainage system 
consists of ditches instead of pipes and the extrapolation of the regression 
equations to Meadowvale or any other area is not justified. The results 
do, however, illustrate the range and variability of storm water pollution. 

Secondly, comparing the SWMM and STORM predictions for 197*+, it 
seems that BOD and SS yields are markedly different (SWMM BOD = 9,400 
lb and STORM BOD = 5,600 lb; SWMM SS = 117,000 lb and STORM SS = 37,300 lb). 
These values have been adjusted for difference in runoff (i.e., SWMM 
yield increased by 1.4 to include missing rainfall and STORM yields 
decreased by 1.1 to account for the suspected overestimate of total runoff). 



TABLE 13- WATER QUALITY DATA FOR CREDIT RIVER. Taken at 
Station C-*t.9 (Highway No. 5), from Ontario 
Water Resources Commission "Water Quality Data" 



DATE 


B0D 5 


SS 


TOT N 


TOT P , 




mg/1 


mg/1 


mg/1 as N 


mg/1 as P 


30/06/65 


1.7 


17 


0.66 


0.60 


14/07/65 


1.5 


6 


0.72 


0.50 


17/08/65 


1.7 


13 


0.47 


0.60 


31/08/65 


1.6 


5 


0.27 


0.70 


6/06/66 


3.6 


15 


1.28 


0.84 


6/07/66 


3.2 


6 


0.43 


1.60 


5/10/66 


1.0 


3 


0.66 


0.20 


18/10/66 


2.4 


9 


1.05 


0.20 


13/06/67 


2.7 


25 


1.63 


0.07 


27/06/67 


3.4 


64 


1.70 


0.18 


31/07/67 


0.9 


4 


0.67 


0.16 


15/08/67 


1.3 


9 


0.71 


0.26 


12/06/70 


5.0 


15 


1.22 


0.24 


30/06/70 


1.4 


10 


0.97 


0.22 


16/07/70 


2.0 


50 


1.57 


0.24 


6/08/70 


1.0 


5 


1.53 


0.24 


10/09/70 


1.8 


15 


2.20 


0.39 


27/05/71 


5.5 


50 


1.47 


0.16 


29/06/71 


2.0 


15 


2.09 


0.24 


14/07/71 


1.2 


45 


1.10 


0.20 


12/08/71 


3.0 


30 


0.99 


0.16 


8/09/71 


1.6 


10 


1.07 


0.21 
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TABLE \k, COMPARISON OF PREDICTED YIELDS USING SWMM, STORM AND 

REGRESSION EQUATIONS (Assumes Meadowvale Fully Urbanized) 



POLLUTANT 


POLLUTANT YIELD LBS (Kg) 


SWMM 1 
PREDICTION 


STORM 2 
PREDICTION 


REGRESSION 3 
EQUATIONS 


BOD 


9400 
(4270) 


5600 
(2540) 


- 


SS 


117000 
(53153) 


37300 
(16950) 


354000 
(160820) 


TOT N AS N 


- 


1500 
(680) 


170 4 
(77) 


I TOT P AS P0 4 




142 
(65) 


480 
(218) 



1 - SWMM PREDICTIONS CORRECTED FOR MISSING RAINFALL 

2 - STORM PREDICTIONS ADJUSTED TO ACCOUNT FOR SUSPECTED OVERESTIMATE 

OF RUNOFF 

3 - REGRESSION EQUATIONS OF COLSTON AND TAFURI 0975) 

4 - ONLY KJELDAHL 



This is a strange result since SWMM and STORM use similar algor- 
ithms for modelling runoff quality. The reason for the difference is almost 
certainly that the "number of dry days" parameter assumed in each of 
the SWMM simulations was too high, giving a conservatively large estimate 
of pollutant buildup. 

The wash-off of pollutants is also affected by the rate of runoff 
as illustrated by the following numerical example. 

For suspended solids prediction, both models use the equation: 

P 0FF = A X P 0N X (1 " e '™ J) 
where: A = availability factor = (0.057 + l.^R ' ), 

p ™ = suspended solids on catchment (assume = 10 lb), 
ON 

P = suspended solids washed-off (the yield), 

K = decay coefficient = A. 6, 

R = runoff (in/hr) , 

T = time (hours) . 

Assuming a buildup of 10 lb of SS on the catchment and an extremely 
simple runoff pattern of 0.5 in/hr for 30 minutes, one obtains: 

For SWMM, R = 0.5 in/hr T = 0.5 hr 

SS 0FF = (0.057 + W x 0.5 K1 ) x 10 x (1 - e^ 6 X °- 5 X °- 5 ) 
= k. 8 lb. 

For STORM, R = 0.25 in/hr T m 1 hr 

SS QFF = (0.057 * M x 0.25 1 ' 1 ) x 10 x (1 V 4 - 6 X °- 25 X K °) 
= 2.5 lb. 
While for BOD predictions, the models use the equation: 
P 0FF =P 0N <' -e _KRT ) + 0.1 x SS^ 

where: P nll = BOD on catchment (assume 0.5 lb), and 
ON 

P = BOD washed-off (the yield), 

substituting applicable values into the BOD equation gives: 

For SWMM, B0D QFF = 0.5(1 - e^' 6 X °' 5 X °' 5 ) + 0.1 x k.B = 0.82 

For STORM, B0D ncr = 0.5(1 - e"^ 6 X " 25 X K °) + 0.1 x 2.5 = 0.59- 
OFF 



Thus, it is seen that STORM underestimates suspended solids 
yields because it underestimates peak flows. Soluble pollutants such as 
BOD and N are not as sensitive to the runoff rate. The STORM results 
discussed previously are acceptable as first-order estimates and the 
comparative analysis is relevent. 

7.3 Meadowvale Reservoir Water Quality 

The success and further application of the scheme at Meadowvale 
will largely depend on the quality of the reservoir water. If continual 
problems of odours, aquatic weed growth, algae blooms, etc., are 
experienced, the developer will not repeat this design strategy. Two 
factors determine the quality of the water, i.e. the characteristics of 
the catchment and the performance of the settling basin. 

The system was simulated using the University of Toronto 

version of STORM. The settling basin has a predicted efficiency 

(Table 15), in terms of suspended solids removal, greater than 85% 

for $6% of the time. This very high efficiency is to be expected 

since with an area of 39,000 sq ft overflow rates are very low. The 

maximum flow rate in the 14-year period simulated was only 209 cfs 

2 
(3^10 U.S. gal/ft /d). The basin wab designed using Stokes' equation 

for the settling velocity, "to settle out all particles larger than 0.07 

mm (fine sand) during a discharge of 407 cfs" (Project Planning Associates, 

Ltd., 1975). 

Pollutant inputs to the reservoir were computed using equations 

described before, and the average annual statistics have been listed in 

Table 11. The model does a continuous simulation of the quality of the 

dam water and provides a frequency analysis of the concentration of 

several pollutants for each month. These analyses for the months of 

June, July and August are shown in Table 16. The following observations 

can be made: 

1) Concentrations of SS, SETL S and BOD are low. This reflects 
the assumptions of the model, i.e., a decay of pollutants to 
specified minimum values between storms. It is possible 
that a resuspension of solids could occur during turbulent 
conditions, while a higher BOD may be generated by the 
ecosystem. 
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TABLE 15- EFFICIENCY REPORT FOR SETTLING BASIN (Based on 
Hypothetical Performance Curve) 



EFFICIENCY « 
SS REMOVAL 


CORRESPONDING OVERFLOW RATE 
(U.S. Gal/Ft 2 /D) (m/D) 


PERCENT i 
OF TIME L.T. 


5 


3315 1010 


0.00 


10 


2710 826 


0.00 


15 


2350 716 


0.00 


20 


2096 669 


0.05 


25 


1920 585 


0.11 


30 


1745 531 


o.n ; 


35 


1600 487 


0.16 


40 


1490 454 


0.21 


45 


1385 422 


0.21 


50 


1298 395 


0.24 


55 


1210 368 


0.29 


60 


1131 344 


0.34 ; 


65 


1061 323 


0.58 


70 


1000 304 


0.71 


75 


860 262 


0.92 


80 


730 225 


1.77 


85 


500 152 


3.90 


90 


o 


100.00 



TOTAL NUMBER OF HOURS SIMULATED = 122554 

HOURS OF STILLING BASIN USEAGE * 3796 

TOTAL VOLUME OF SOLIDS REMOVED = 5997 Ft 3 = 170 rrf 
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TABLE !6. MONTHLY FREQUENCY ANALYSIS OF POLLUTANT CONCENTRATIONS IN RESERVOIR, STORM PREDICTIONS 





MONTH NO 


= 5 










RESERVOIR QUALITY ANALYSIS 








S.S. 




SET.S 


BOD. 


NIT. 


P04. 


D. 0. 






CONC 


% TIME 


CONC 


% TIME 


CONC 


% TIME 


CONC 


% TIME 


CONC 


% TIME 








MG/L 


LT. 


MG/L 


LT. 


MG/L 


LT. 


MG/L 


LT. 


MG/L 


LT. 








5.00 


1.937 


2.00 


100.000 


2.00 


100.000 


0.20 


0.0 


0.10 


0.0 


1.00 


0.0 




10.00 


100.000 


5.00 


100.000 


4.00 


100.000 


0.40 


0.0 


0.15 


6.718 


2.00 


0.0 




20.00 


100.000 


10.00 


100.000 


6.00 


100.000 


0.60 


0.0 


0.20 


6.718 


3.00 


0.0 




30.00 


100.000 


15.00 


100.000 


8.00 


100.000 


1.00 


0.586 


0.25 


6.718 


4.00 


0.0 




50.00 


100.000 


20.00 


100.000 


10.00 


100.000 


2.00 


6.718 


0.30 


9.931 


5.00 


0.0 




75.00 


100.000 


25.00 


100.000 


12.00 


100.000 


3.00 


31.210 


0.35 


33.336 


6.00 


0.0 




100.00 


100.000 


30.00 


100.000 


14.00 


100.000 


4.00 


57.120 


0.40 


42.493 


7.00 


0.0 




150.00 


100.000 


35.00 


100.000 


16.00 


100.000 


5.00 


100.000 


0.45 


57.120 


8.00 


0.0 




200.00 


100.000 


40.00 


100.000 


18.00 


100.000 


6.00 


100.000 


0.50 


85.807 


9.00 


0.803 




900.00 


100.000 


75.00 


100.000 


25.00 


100.000 


10.00 


100.000 


0.60 


100.000 


16.00 


100.000 




MONTH NO 


■ 6 










RESERVOIR QUALITY ANALYSIS 








S.S. 




SET 


.s 


BOD. 


NIT 


a 


P04. 


D. 






CONC 


% TIME 


CONC 


% TIME 


CONC 


% TIME 


CONC 


% TIME 


CONC 


% TIME 








MG/L 


LT. 


MG/L 


LT. 


MG/L 


LT. 


MG/L 


LT. 


MG/L 


LT. 








5.00 


1.840 


2.00 


100.000 


2.00 


99.864 


0.20 


0.0 


0.10 


0.0 


1.00 


0.0 




10. UO 


99.981 


5.00 


100.000 


4.00 


100.000 


0.40 


0.0 


0.15 


5.713 


2.00 


0.0 




20.00 


100.000 


10.00 


100.000 


6.00 


100.000 


0.60 


0.0 


0.20 


8.191 


3.00 


0.0 




30. OC 


100.000 


15.00 


100.000 


8.00 


100.000 


1.00 


0.0 


0.25 


8.191 


4.00 


0.0 




50.00 


100.000 


20.00 


100.000 


10.00 


100.000 


2.00 


8.191 


0.30 


15.666 


5.00 


0.0 




75.00 


100.000 


25.00 


100.000 


12.00 


100.000 


3.00 


22.512 


0.35 


23.606 


6.00 


0.0 




100.00 


100.000 


30.00 


100.000 


14.00 


100.000 


4.00 


58.801 


0.40 


38.033 


7.00 


0.0 




150.00 


100.000 


35.00 


100.000 


16.00 


100.000 


5.00 


100.000 


0.45 


59.072 


8.00 


0.0 




200.00. 


100.000 


40.00 


100.000 


18.00 


100.000 


6.00 


100.000 


0.50 


78.495 


9.00 


17.942 




900.00 


100.000 


75.00 


100.000 


25.00 


100.000 


10.00 


100.000 


0.60 


100.000 


16.00 


100.000 




MONTH NO 


= 7 










RESERVOIR QUALITY ANALYSIS 








S.S. 




SET 


.s 


BOD 


t 


NIT. 


P04. 


D. 0. 






CONC 


% TIME 


CONC 


% TIME 


CONC 


% TIME 


CONC 


% TIME 


CONC 


% TIME 








MG/L 


LT. 


MG/L 


LT. 


MG/L 


LT. 


MG/L 


LT. 


MG/L 


LT. 








5.00 


1.741 


i-J 


99.991 


2.00 


99.944 


0.20 


0.0 


0.10 


0.0 


1.00 


0.0 




10.00 


99.963 


100.000 


4.00 


100.000 


0.40 


0.0 


0.15 


0.0 


2.00 


0.0 




20.00 


100.000 


10.00 


100.000 


6.00 


100.000 


0.60 


0.0 


0.20 


7.131 


3.00 


0.0 




30.00 


100.000 


15.00 


100.000 


8.00 


100.000 


1.00 


0.0 


0.25 


7.131 


4.00 


0.0 




50.00 


100.000 


20.00 


100.000 


10.00 


100.000 


2.00 


7.131 


0.30 


13.945 


5.00 


0.0 




75.00 


100.000 


25.00 


100.000 


12.00 


100.000 


3.00 


26.382 


0.35 


27.090 


6.00 


0.0 




100.00 


100.000 


30.00 


100.000 


14.00 


100.000 


4.00 


65.435 


0.40 


34.742 


7.00 


0.0 




150.00 


100.000 


35.00 


100.000 


16.00 


100.000 


5.00 


100.000 


0.45 


66.375 


8.00 


0.0 




..|M°.. 


..188:888.. 


J8JL 


.ML 


Ji;8L 


100.000 

_1Q0,00D- 


6.00 


100.000 
.10Q.Q0Q-. 


0.50 
_.Q.6Q_-. 


88.801 
.100.QQ0- 


9.00 
._16.QQ__ 


100.000 
.100.000-- 






2) Concentrations of nutrients (assuming no sources or sinks) 
are significantly larger than the threshold values (0.3 mg/1 
and 0.01 mg/1 for inorganic nitrogen, N, and inorganic 
phosphorus, P, respectively) reported by Sawyer (1952) for 
the occurrence of algae blooms. Even if some nutrients 

are removed (e.g., by the growth of aquatic plants or the 
formation of complex precipitates), conditions suitable 
for algae growth will likely prevail. There are various 
alternatives for algae removal; one which appears suited 
to Meadowvale is the use of dilution water (Rohlich, 
1969). A flow of I cfs during the months of June, July 
and August would reduce the average phosphorus concentration 
from 0.^2 mg/1 to 0.16 mg/1. The fairly fast turnover of 
water (once per annum) and the shape of the reservoir 
(conducive to plug flow) will tend to wash algae out 
and the problem will not necessarily occur. 

3) Dissolved oxygen levels are close to saturation. As out- 
lined, the decay of vegetation or algae could place a demand 
on the oxygen resource, however, depletion below 3 to k mg/1 
is unlikely and the reservoir should provide an adequate 
fish habitat. 

k) Should anaerobic conditions develop in the settling basin, 
there is potential for a shock loading on the reservoir. 

7. h Sediment Yield and Erosion 

Urbanization generally leads to an increased sediment yield from 
a catchment and to accelerated erosion in the plain of the river or 
stream which receives the urban runoff. Most of the sediment yield is 
due to smaller storms. Glymph ( 1967) studied the records of 72 water- 
sheds in the United States ranging in size from 100 to 100,000 acres 
and concluded that storms with a return period of less than one year 
contributed between 3^ and 92 percent of the aggregate sediment yield. 
As discussed before, the Meadowvale system can be expected to remove 
about 90% of the incoming sediment and in addition will reduce peak 
flows, which are known to be the major factor in river bank erosion. 
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For urban areas the major cost of erosion is not the loss of soil, but 
rather the cost of removing silt and debris from estuaries, reservoirs, 
and stream beds (e.g., the sediment deposited by the Don River in the 
Keating Channel of Toronto Harbour necessitates annual dredging of some 
100,000 cubic yards) . 

Construction activities during the transition from rural to 
urbanized conditions are a major source of sediment. Yields as high as 
140,000 tons/square mile have been reported (Dawdy, 1 967) for small 
construction sites. The SWMM model has the capability of modelling erosion 
using the Universal Soil Loss Equation. An attempt was made to determine 
the order of increase in sediment that could be expected in Meadowvale 
during the development of land to the west of Winston Churchill Boulevard. 
These results are given in Table 17 and indicate a four-fold increase in 
sediment yield (i.e., from 60 tons/year to 240 tons/year) during the 
construction phase. This corresponds to a yield of between 1 60 to 
640 tons/square mile/year. Since this part of the catchment will not be 
developed gnti 1 after the dam is complete, the "computer guesses" above 
could easily be verified by the Ministry of the Environment in their 
program to monitor Meadowvale. It is unfortunate that the Meadowvale 
dam was not built prior to the construction of all the houses, roads and 
storm water sewers. Ideally, if a siltation basin storage system is 
planned, it should be constructed first so that maximum benefits may be 
obta i ned. 
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TABLE 17. PREDICTED INCREASE IN SEDIMENT YIELD DUE TO CONSTRUCTION (SWMM) 



STORM OF MARCH 4 

TUTAL lU'„Pt Nl)l T) SCC IRS (INCLUDING 



.. f-.JH I H;',li;M riiANbFI.ur.TD to THf- SFWls rni)N tNLFT ARFAS DURING 7« TIMf 



INL " (lij'I IN FLOlK ','HL CHOPPING CUNTRUI. 

no. »!i»(«cmr.i it'ir.ihirn iyi'Kki facto' FAcmi- 



fl'l]:. I UN 
s.s.i mgi 



TOTAL 
S.S.I MG> 



1 

4 
S 
G 

r 
ft 
9 
II 

1 1 



0.50 

.HO ' 

O.fiO 

ib.L.0 

■..70 
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I 3. 00 

?5.00 

15.00 

ss.ee 

1 .00 
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~~ /CO. 00 
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APPENDIX 

SUBROUTINE "MIKE" - 

PROGRAM SOURCE AND MODIFICATIONS TO STORM 



ON 



is 

10 

21 
11 

13 
34 
15 
16 
17 
It 
1» 
ID 
11 
33 
33 
34 
35 
36 
37 
J8 
,3» 



SUBROUTINE 

COMMON 



MI KE 

I PAIN (200, 24),JDATE(2C0) . JDAY(20C ) . 
KRAIN(200.24>.KDATEC2C0).KDAY(200> 



COMMON /REEL/ 



1 
2 



AGE, AGE 2, AGE 3. AGF4i AGE 5 , A GE 6 .AGE 7 • ACE 8 . 
AMXOLD.CC IAT. CIMP.CN.CPEPV. DD 1 . HEP , DE P K , DEPRF. 5» 
DEPRN .DEPRS , DUR ,E ,EXPTE .EXPTN ,HUND. OLD . FRCP. 



-~5~ 
4 
5 



P R f .P N , '3 ft INR1 ,frA I NX «PF TF ,R F N . P F U . PT OlT S , F UN t 

SUNOP ,SACFT,SAGE.5AGE2,SAGE3,SAGF4,5AGE5.SAGE6. 
SAGE 7 .SAGlP . S«f, , SPUN , T ARE A ,TCFS , T MGD ,U R C, 
Xi\GE, XA Gt 2.X AGE 3.X AGE 4, XAGE5, XAGE^.XAGE 7, X AGE 8. 
Xi-'UN'i YP A IM .Tf-'UNi 7EX 



COMMON /DIMN/ 



CAPRI tO »S0 ) , CP ( 5> , QD( 5) ,DELTP( 5, 5) ,DELTPN( 5) . 
FIMP(S) .Ft- ACTN(5,5>,P(5.5).PN(5) . POL (5 ) ,°OLEX{ 5) . 



P 
3 



P'JLLKT ( 5) 
REC.VRTt 12 



PCLNRTI5") iPDLOVR( c ) .RECVNt 121 



> ,5SUMP0< 5) , STLEN(5 > ,SUMPOL< 51 . 
TrATEI<(20) ,XP0LEX(5) ,XP0L0V(5 ) .XSUMPC(5 > 



CUMML'N /I TGRV LAP, LSI . tX ,hV(_bS, 1A. I <3 , ICC, ID, I DATE, IENH, IL. IN 



A 

n 
c 



ISCH.IT , 



F 
F 
G 



!"«? . IPiCK, I PR INT. I PI . IP2 .IP . I RUN, I ST APT , 
J, JHn?,JMX,KK,KMAX2,KMX,L.LAST,LPATF , 
LEXT, L MUNO ,L I , LL.LCSS , LSDF , M , MAS E . M ASS , 

TTA XLT N, M; 



,iviK«^t .i«ST, VXC. 

MXTIM.NAM.NANTEC.NCtNCAP.NDAYPC.NDAYl.NrrPV, 
NHR,NTNC,NMAX,NN,NNRAIN,NCEXTRtNCRAIN, 
N'^PUN .MQnUNN .NOSTCR .NPAGE , NSUIVR . RA If>' , SO AT£ , 



MXSTOP, 



~FT 
I 
J 
K 

T7 



'j l JAV . VVJM . SEX . bLXLFb. SHK, 5W?T, 5 M T S TP . 

^NC«TC,SRAIN,STCAP.STNSTC.STC D .STCR, 
STGRQ .ST'OPG-j STRAIN, STRE AT.STSTO.TCAP, 
" r 'JST'?r;,TCP,TPAIN.TREAT,TSTOR.XCST,XPU D . 
X" A I N . X TNS TO . X TOPT 



5500 
5501 
5503 



5 5 04 
5505 
5506 
55C7 



5 508 
5509 
5510 
551 1 



5* 12 
5513 
5514 
5515 



55 15 
551 7 
551 * 
5519 



5 520 
5521 
5522 
5523 



5525 

5526 
5527 



5 5? 8 
5529 
5530 
c 531 



COMMON /IOIM/ 



1 

~2~ 



DSTQf; (^00 ) . ICNA'JtE < 5 > . IPOLUT ( 20 > .LPAGE (5 >. 
LSniSTC 100 ) .LTDIST( ICOO) , NAME (8) . NAMEWS( 4 ) , 
NCLhAM t 5 ) , NP ( 12 ) .N TIT LE t ? C > . TP A T E t 2 C ) . Y« UMPf ( 



T 



5533 
5 5 34 
5535 



5^» f i 
5537 
5 5 38 
5539 



COMMON/ U0FT/PV3L ,AR»1 . WTP2.F- VOL . NHTS ,HAJ(5),ARJ(5),0C1,S1, 
U WOI^.CONC f; ) ,PPJT2( 51 ), DONE, STORP( 51 ) ,CFLC*( 51) ,°REV, 

SU RE . S C r NCtP . TOT , SI W f 6.1 .12 ) . J P>PEV,LX ,-D&t t 12 » . 
TEMPI 12 ) , 3EF F ( 10) , MCOUNT, SOLIDS ,M£FP 
DECLARED INTEGER VARIABLES IN CDMMON 
INTEGER CAP,CST,DSTOR,DUO f tx,fcXCE5.RAIN,SDATE,SDAY,SCUR, 






5-540- 
554 1 
■^542 
554 3 




C 554* 

COMMON /SOIL/ A-3F A, AP EAN, TER , IFPDMX , IPRTS , I3.L INE.MCF ,m I 30.MXLG, 5549 

1 NONURB. NSTOR .REC.RMI ,SMEC , XMH, 5550 

2 DDL (6 ), EI( 6 ) , I EPP < 20 ) .LNDUSF « f- .2 ) iPPCM (6 ) ,WTQ<6) 5551 
C 55^2 

INTEGER URE 555"} 

t/« REAL OVFL, S202 •S2.Q02 .POL IN( 6) 5554 

/T 150 0C FORMAT! » I • . 1 J2T. ' *' J ,/ ,4*X , » WLSEPVLIW hbi'UVl — HJR bU bt I M) E Vt N I S ' . / bbbb 

A.SX.'EVNT'.TX.'DATc'./X, 'HH ' ,6X, 'f)I N' »7X,' STOP AGE" . 6X ,• QCUT * »6X . 'SS556 
B.S.* .5X , 'SETL.S. • .5X, • PCD* . 7X . • TGT N * , 7 X . • P • , 7 X . *0 . 1 . • ) 5557 

3 15100 FORMAT! 1 I X ,« YR MQ D Y» , 1 3 X , • C FS' . &X , * F T* * J • , 7X . » CF5 • , 7 X . • MG/L * , 6 X . • 5558 

— CMG/L'.TX, »MG/L' . &X. 'MG/l • .f X , 'H«/L' 16X1 •H(,/L« ,/ . I.«7t * »* H ""5559 

15200 FORMAT(I6.4X.3I?,5X,IP,5X,F5.1,5X,F9.0.5X,F5.1,3(Sx,r' : .C).l(5X,F-5.55^0 
D2) > 556t 

6 MPRNT = Q 5562 

, I NSTDP = NSTDR+ 1 5"5"F3 

IF{ IPOLMX.FQ. 0) GO TO 11115 55631GC 

DO 11110 K=1.IPHLMX 5564 

IF ( IPOLUT( Kl.tQ.INSTOP.CR. IPOLUT(K ) .EQ.( INSTC&- 1 ) ) MPRNT=1 5565 

"TTI10 COrTTTNTJE - ^565^ 

11115 CONTINUE 5566100 

IFfOONE.GT.l ) GO TO 11160 5567 

C FIRST POUTING OF FLOWS 8568 

C CALCULATE — 5 1 AGE STORAGE DTSl_HAHGb AT — DUHtPLN'l — HL A V \, Sbfcy 

LX=0 55/0 

<=1 5571 

SOL10S=0. 5572 

HEFF^TJ ^.^m 

PPEV=0 557* 

JPREV=24 5575 

MCOUNT=0 =5 5 76 

ROUT?! 1 ) = 0. ^b/; 

STORR1 1 )=0 . 5 C 7P 

HEAO=0. 5579 

2' OFL3W( I )=0. 5580 

„ — "S2=TJ. 55^fr 

001=C. 55811CC 

" 00 1 1 150 N=l, 50 55B2 

i* HEAD=HEAD+0 .1 5 r ei 

IUHfcAU.tiT.IHAJ^.'HlM)) Tji: !L 1116 55^U 

I F( HEAD .GT . (HAJ (K ¥\ ) ) } K=K*1 ^585 

STORPIN + 1 >=STCRP(N»+0,1*( ( t-EAD-FAJ(K )-C .05>*{s(~J<K + l |-AOj(< ) ) / ( HA J55P6 

.(KH)-HAJ(K))tAKJ(K|) 5587 

~tftwtp27 — n 12c, i li ic. i i l^o 5^8^ 

C OGEE SPILLWAY 5590 

1 1 1 2C OFLOWtNfl }=?. ?*wniM*HtAO**l. 5*3600. 5591 

GO TO 11145 55^2 

C MIJHN1NG GLLMV "^93 

1 1 I 3C QFLOW (N«-l) = (4.3 5-.l Q* HEAD /WD IM)*2.*3. 1416* WPT.M* HE AP**1. 5*3600. b^ 94 

" GO TO 1 1 145 5 5C4 

53 C CONTROLLED CONSTANT DISCHARGE 5595 

1 I I 40 - TR. O W" TTJ> 1 )=36C0. *WD I M =? 5<?6 

11145 CONTINUE 55<?7 

3S ROUT2I N + 1 > =OFLOW( Nt- 1) +STOPP(N+l ) 5598 

" 11150 CONTINUE 55Q9 

— r yuy^tNUbn sLLiub (aieggptes if c o 

SCONCC 1,1) =5. 5^01 

M SCONCI 1 ,2 )=10 * 5602 

39 SCONCC 1 .3 >=20. 5feC3 

sconci i r4-)=ju. --wtr*— 



10 



u 



3? 



SCONCf 1 ,5 >=SO. 56C5 

SCONC ( 1 .6 ) = 75, 5606 

SCONCf 1 , 7)=10C. 5607 

SCONCf 1 .8) =1 50 . 5608 

SCCNC ( li<i>=20r, 560O 

SCONCf 1 , 10 )=^C00 . 5610 



1/6 


c 


5ETTLEABLP SCLfOG CATEGORIES 


^.61 1 


A 




SCONC f 2. 1) = 5. "'" "'" 


^fc It 






SCONC( 2,2) =1 0. 


561 3 


i 




SCDNC (2,3 )-15 . 


5614 


3 




SCONCf 2 ,4 )=20. 


5615 




SCONCf 2,5 ) = 25. 


5bl 6 






SCDNCt2.6 J=^0 i 


561 7 


1 




SCONCf 2,7 ) = 3^. 


5618 


t. 




SCONCf 2.8 > = 4Q, 


5613 




5C0NC(2.9 t=Ab. 


b620 






SCONCf 2,10) =75. 


5621 


A 


c 


POD CATEGCRICS 


5622 


9 




SCONC (3 . 1 ) =5. 


5623 




SCONCf 3.2)^1 C. 


5624 






SCONCf 3,3 )= 15. 


5625 


1 1 




SCONCO ,4 )=2C . 


5626 


13 




SCONC( 3, 5)=?^. 


5627 




SCDNC I 3,6 )= -, C . 


^^>cLU 






SCONCf 3. 71-4C . 


5628 


14 




SCONC f 3, 8 1 = SC . 


5 62 <=> 


IS 




SCONCf 3,9 )-? K . 


5610 






sc owe c 3 , ro \ =1 e& . 


5631 




c 


NITPDGEN C ATFC.~- ; ] r S 


5633 


1 ( 




SCONCf 4. 1 ) = C «2 


56 34 


1. 




SCONC f 4,2) =0. 4 


5^35 






5C0NC(4, I!=0,^ 


56.4 6 






SCONC (4 ,4 ) =0 , c> 


5637 


10 




SCONC ( 4,5>-l . C 


5638 


II 




SCONCf 4,6 >=2 »C 


5639 




SC0NCC4 ,7) =4 .C 


5640 






SCONCf 4.8)=8. 


5641 






SCONCf 4,9 ) = 10 . 


5642 


34 




SCONCf 4,10 )=?C .0 


5643 


l?l 


c 


PHOSPHOROUS L"JM hi-i'M ltflb 


5644 


1 




SCONCf 5, 1 »*@»G2 


5645 


u 




SCONCf 5 .2) =0 .04 


5646 


37 




SCONCf 5,3) = 0. C6 


5647 






SCONCf 5,4 |*0 .OR 


56T8 






SCONCf 5,5) = C. 1 


564*9 


2? 




SCONCf 5,6) =C. ? 


5650 


30 




SC0NC(S,7 > = 0.4 


5651 






SCUNCf biH)=0,B 


3e32 






SCONrf 5,9>=1 .C 


5653 


1? 




SCONCf 5 ,10 >=5 .0 


5654 


13 


c 


DISSOLVED OXYGEN CAT^GCRItS 


5654 






SCONCf 6, 1 > = 1 • 


5655 






SCGNCf6 ,2 )=2. 


5656 


3S 




SCONCf 6,3>=3. 


5657 


36 




SCONC(6,4 )=4. 


5658 






SLUNL I b ,5 )=b . 


5 6 5 9 






SCONC (6.61=6. 


5 66 


38 




SCONCf 6, 7) = 7. 


5661 


l" 




SCONC(6 ,8 )=P . 


5662 


|» 




" SCONC I b ,9 ) -9. 


Son J 



SC0NCI6. 10 > = 1 f • . 
0021 1 ) =14.6 
D0 2( 2 1 = 14.6 
D02f 3 > = 14 .6 
0021 41=12.4 
0021 5)= 1 0.8 



56 64 
5665 

5666 

566? 
5669 
5 6*9 



i o 




D02 (6 )=9.5 


56 70 




fi 




DG2I 7J=a.l 
DO 2 C B ) = 9 . 2 


5R m 

567? 




J 




D02{ 9 ) = 10 . 1 


5673 




i 




D02( 1 t >=13.3 


56 74 






D02( 121=14.6 


"56T5 








OXDEF=D02I MO)-CONC(ft) 


5*76 




1 


C 


MEAN MONTHLY TEMP 


5677 




t. 




Tc«P( 1 »=4. 


5673 






TEMC>| ?) =4. 


56 ry 








TEMPI 3) = 6. 1 


5680 




a 




TEMPI 4 » =6.4 


5681 




9 




T=MC( 5)=12. 


56P2 






TEMP( 6)=19. 2 


^56R3 








TEMPI 7>=20.7 


5 6 84 




1 I 




TEMPI 8)=20.0 


5 6 85 




11 




TEMPI 9)=15.7 


5686 






TEMPI 10 J=9 .8 


568/ 








TEMPI 1 1 ) =4.0 


5 6 88 




14 




TFMPI 12 »=4 ,0 


5689 




11 


1 11 60 


CONTINUE 


56QC 






LX=LXM 


56 91 






1 1180 


CONTINUE 


5 6 3? 




1/ 




HEFF=HEFF+1 


5693 




18 




RUNOF=RUNOF*PT0CFS 


5 6Q4 




C 


ACCOUNT FOR 1 I Mfc 1NILWVAL HLIWbbN KUNl't-1- S 


5 £ 9S 








IFIKDAYILX ) . GT.PTEV > GO TO 11190 


5^96 




M 




UR£=J- JPREV-1 


5 6 97 




>l 




GC TO 11200 


5698 




I1T9TT 


URfc = 2 3-JPHEV+ (KIJAYI LX J-PREV-I )*24 4- J 


5699 






c 


URE IS THE NO OF HOURS WITH NO RUNOFF 


57C0 




13 


1 1200 


CONTI NUE 


5701 




M 




SUPE = SURE + URE 


570? 






II- lUKh .bU.U ) r,U 'U 11355 


5703 








00 1 1 350 N=l , URr" 


57C4 




It, 




S1=S2 


5 705 




J' 




S202=S 1-OC1/2 . 


5 706 






IF( S2C2.yi .0 . 1 GC TO 1121T 


^70T 








S2=S1-PECVRT( *0) /28P.*AR J I 1) 


S708 




20 




Q01=0 . 


5709 




JO 




002=0. 


5710 




c 


Nl) OUII-LUW WE SEW VOIR NOT FULL 


57 I 1 








GO TO 1 1305 


571? 




VI 


1 12 10 


CONTINUE 


5713 




yj 




DO 1 1250 K=l, f-0 


5 714 




1 


5T1 5 






1 12 50 


CONTINUE 


5716 




Ji 


1 1260 


S2 = ST0RRIK>-RFCVRT| MO 1/288. *AR J( 1 ) 


571 7 




30 




0O2 = Cf r LOWl K > 


5718 








aoi=ou2 


5 719 






113 


CONTINUE 


5720 




18 


1 1305 


CONTINUE 


5721 




19 


C 


UPDATE WATER DUALITY U5ING DEFINED CATEGORIES FOR *=P C C ANALYSIS 




M 


C 


KZ1 IS UEDXYGENAT ION COEF KZ2 IS REACRATION CPEF 


3 l S. J 





I 6 



KZ1=C . 1 *1 . 04 7** ( TrS'Pf MCt -2C. ) 5724 

K72=0.0 75»l.C»'***( T iMP(f"C)-2C.) 5 72 5 
FAT l=K71/( K71 -K /.? ) 5726 

0XOEF = FAC 1 *C3fcC < 31* H 0.**<-K71*N/24. ) -1 0. ** (-KZ2*M/24 . ) ) taXDEF*l .S7?7 
, **(-K72*TDAYS > S7?8 

CONC ( ft ) = OQ2 ( WC > -Txr)[_r 5729 

CONC ( 3)-C0NC( 3) * 10. *v ( K7 1 *N/ 24. ) 5730 
DO miO K = l , 10 577T 



TF (CONC ( 3) .LT.5CCNC (1 .K)> GC T0 11315 57"<2 

11310 CONTINUE 5733 

11315 CONTINUE 5 734 

SIM( 3, K, MC T^S I M{ 3 »K,MBJ +.3 5 735 

DO II 330 K=l, 1C S7^6 

IF(CCNCto) .LT.SrpNC (6«K» ) GC TC 11335 5737 

11330 CONTINUE 5738 

11315 CONTINUF WT&9. 

SI V( 6, K,MC) =5 TM(ft ,K ,MC ) + l 5740 

DO 1 1 336 K=l , 1 C 5741 

IF (CCNCI 1 ) .LT .SCONC C 1 iK > ) 5 I M ( 1 . K , MO ) = S I M ( 1 . K , MO ) + 1 5742 

IF I CUNC ( 2 I •ITTVSC'JNiL » «£ » K i I — STWTSTTK .I^TOT^S TT»TT2"i K rMTTTrl 5T4 3 

11336 CONTINUE 5744 

" CCN<" ( 1 )=CCNC ( 1 ) *-xr>»* (-1 ,,?C4 ) 5745 

CONC ( 2) =CCNC ( 2) «CXi»«* ( -1 . 2C4 > 5746 

TTICCNCI 1 ) .1 r .'-. If -"M ( 1 ) -b. 574/ 

IF (CCNC( 2) .L T .2 . > C£TNC(2J=f« 5 74 8 

14 IF(MPRNT.rO.O) G IT C 11340 5 7 4Q 

WCCUNT =MCCUN T « 1 575C 

I FT MC UK'T . L T."^T 1~ T^T" 'T -- ri~T3R " 575 I 

«o I T »- ( 16, 1 50CC ) 5752 

W P I T F ( 1 6 , 1 5 1 :. ) 5753 

MCOUNT = l 5 754 

1133P CONTINUF 5/55 

WPTTE( 16 , I520C > MSTOR ,kvc , km n , KOY , J.PUNCF, 5 2,002. ( 03NC( K|_> , KL = 1 .6)5756 
1 1340 CONT INUF 57=^7 

11350 CONTINUE 5758 

JPREV = .T 5 7"^P 

PRE V=KOAY ( I. X ) 576C 

S IM( 4 ,K 14 , MO > =S! «< 4 ,K 14. MQ)+-URF 57M 

a< SXWC 5»ttl5»«0»*S£'M|-S tKl^i MO> +URE 5762 

C N I' OW HOU TE — FTT^ — p-"^ — FT-~TT7~TT FTJNTJFT ? 763 

11355 CONTINUF 5 764 

POWFC =.000285 5765 

37 HEFP = HEF tr +■ 1 5766 

DVFC«l?UNO*-/AHMl *64fe .*1 /» Br7WT 



■i 
i 

« 

a 

9 

10 



IF( OVFL.GT, 1 CCC . ) SC TO 114CC 5768 

" FFF=. <J- ,0C 07*PC *K P*C VFL* *2. 5769 

jo GO TO 11410 5770 

ri4oe ErF = c. 7 %"x p ynt ' ?'iti-t-' » 2.«i louc .-cv ^ l ? ) ^tt 

1 1410 CONT INUF 5772 

33 DO 1 1 440 N»| , 10 5773 

STFP=N*10. 5774 

tf~ r ?TFP".i3"r.rFFT (sc to ~tt 445 - ~ trr 75 

I 1440 CONT 1 NUE 5776 

35 11445 CONTINUF 5 777 

K17=N 5778 

Sb>t- I Kl / } =bt-l-V ( R 1 7 ) +■ 1 ^779 

POLLUTANT INPUTS IN POUNDS 5780 

31 P0LIN< 1 )= I 1 .-EiFF )*=>OLLRT( I ) 57B1 

P0LIN(2)=(1.-FFF)*PlLLPT(?) 57 82 

P0LTNf3)=( r.-T.5*EFF-) HP OLLR TI 3 I -»T<«I 



," 



POL IN (4 > = ( 1 ,-f, . TwFFT ) *PQLLPT {4 ) 
POLIN( 5»=( 1 ,-C.2*£FF) *POLLPT ( 5) 
POL IN (6 ) = 2 ,0?3*9UNDF 
SOLI DS=SCL [DS+?Pr*oCLLRT(l )/S0, 
S20 2= 36 00. *RUNOF +-51 -0 1 
IF(S2C2.LT .0 . > GO TO 11456 



57^)4 
57P5 
5786 
57P7 
5788 
5789 



« 




DO 1 1 45C K = l , 50 


579C 




f, 




~TF'OnUT2tK ) .GT.S202) GO TO 114K3 


5 7^1 






1 1450 


CONTINUE 


5792 




t 




GO TO 11460 


5793 




3 

.1 


C 


PESFRVOIR STILL NOT FULL 


57R4 




"IT 455 


52 = 51 «-RUNOFF*3600 . 


^T<?5 








002=0. 


5796 




5 




001=0 . 


57Q7 




6 

7 




GO TO 11475 


5798 




1 1460 


52=5TORR[K J 


5 7Q9 








QO2 = CFL0W (K ) 


5800 




e 




001=00? 


5P01 




10 


C 


POLLUTANT ACCOUNTING FOP THE RUNOFF 


5802 




~I 1 4T5 


CTJNTTNLTF 


SP02 








DO 1 1 480 1=1.6 


5P04 




1 1 




CONC ( I ) = (CONr ( I |*(SUHVCL) *2 8. 32 <-POL I Nf I |*4S35 ( 5?.4 >/2P. 32/52 


580"=; 




ii 




IF(LX .E0.200 ) LX-0 


5 806 




1 14R0 


rUNI INUE 


bye T 








IF( MPRNT.F. O.O 1 GO TO 11490 


5808 








N*COUNT = MCCUNT + l 


58C9 








IF(MCnuNT.LT.bO) GO TO 11485 


5810 






*r itec 16, i sooirr 


sBir 








WRITE( 16, 15100 ) 


5812 








MC0UNT=1 


5813 






I 14 85 


CCNT T NUF 


5814 






WRITE! 16,15200) N5T0R.KYR , KMO.KOY , J.WUNGF, S2.UU2. [ LJ N £ ( KL > . KL = 1 


, fc > 5 H I 5 






1 1490 


CONTINUE 
SI = S2 


5816 
5 81 7 








DO 1 1 540 1=1.6 


581 8 






TO 1 1 530" K D~= 1 . I 


5BT9 








IF(SCONC( I ,KD) .GT.CCNC( I ) 1 GO TO 11535 


5820 






1 1530 


CONTINUE 


5821 






I 1 5 35 


CONTINUE 


^822 






b I M( I ,K!),MUJ=blMl l.KD.MUJ + l 


56^3 








if( i .eo. i j ki i««n 


5624 








IF( I .EO.2 ) K 1 ?-*V> 


5«?5 




ik 




IF( I .EQ.3) Kl 5=KD 


5826 






IF*I«fcO,4J K 1 4 = Kf) 


5°?7 








IF( I .EQ.5I Kl 5=Kn 


5 82 8 




J« 




IF( I .EQ.6) KI 6 = K0 


5 829 




JO 
3 | 


1 1540 


CONTINUE 


5P30 






CUMC 1 1 J = CUNC I 1 > *; Y.P ** (- 1 ,W* ) 


^8^1 








CONC ( 2 ) =CCNC( 2 ) *FXP** (-1 ,2 04 ) 


5 8.32 




11 




IF(CCMC(1 I.LT.S, )CGf C ( 1 ) =5. 


SB? 3 




33 




IF (CCNC (2 > .L T .2. ) CCNC(2» = 2. 


5 6 34 








RETURN 


?rf*~*^ 





END 



5836 



